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Capillary  pressure  vs.  saturation  (Pc(SO)  curves  are  fundamental  to  understanding  liquid  water  transport 
and  flooding  in  PEM  gas  diffusion  layers  (GDLs).  Pc(Sl)  curves  convolute  the  influence  of  GDL  pore  geom¬ 
etry  and  internal  contact  angles  at  the  three-phase  liquid/solid/gas  boundary.  Even  simple  GDL  materials 
are  a  spatially  non-uniform  mixture  of  carbon  fiber  and  binder,  making  a  Gaussian  distribution  of  contact 
angles  likely,  based  on  the  Cassie-Baxter  equation.  For  a  given  Gaussian  contact  angle  distribution  with 
mean  (0Mean)  and  standard  deviation  (or),  a  realistic  Pc(Sl)  curve  can  be  computed  using  a  bundle  of  cap¬ 
illaries  model  and  GDL  pore  size  distribution  data.  As  expected,  computed  Pc(Sl)  curves  show  that  0Mean 
sets  the  overall  hydrophilic  (0Mean<9O°)  or  hydrophobic  (0Mean>9O°)  character  of  the  GDL  (i.e.,  liquid 
saturation  level  at  a  given  capillary  pressure),  and  cr  affects  the  slope  of  the  Pc(Sl)  curve.  The  capillary 
bundle  model  also  can  be  used  with  (0Mean,  cr)  as  unknown  parameters  that  are  best-fit  to  experimen¬ 
tally  acquired  PC(SL )  and  pore  size  distribution  data  to  find  (0Mean,  cr )  values  for  actual  GDL  materials.  To 
test  this,  pore  size  distribution  data  was  acquired  for  Toray  TGP-H-090  along  with  hysteretic  liquid  and 
gas  intrusion  capillary  pressure  curve  data.  High  quality  best-fits  were  found  between  the  model  and 
combined  datasets,  with  GDL  liquid  intrusion  showing  fairly  neutral  internal  surface  wetting  properties 
(#Mean  =  92°  and  a  =  10°)  whereas  gas  intrusion  displayed  a  hydrophilic  character  (0Mean  =  52°  and  cr  =  8°). 
External  liquid  advancing  and  receding  contact  angles  were  also  measured  on  this  same  material  and  they 
also  showed  major  hysteresis.  The  new  methods  described  here  open  the  door  for  better  understanding 
of  the  link  between  GDL  material  processing  and  the  wetting  properties  that  affect  flooding. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Polymer  electrolyte  membrane  fuel  cells  (PEMFCs)  show 
promise  as  power  sources  for  a  range  of  applications  owing  to 
their  high  power  density,  low  pollution,  and  low  noise  [1,2].  How¬ 
ever,  water  management  continues  to  pose  a  considerable  challenge 
[2-6].  When  the  PEMFC  is  operated  at  high  current  densities,  liquid 
water  floods  and  blocks  access  to  the  catalyst  sites  from  the  reac¬ 
tant  gases  [1,4-8].  Part  of  the  solution  to  liquid  water  management 
is  the  design  of  the  porous  gas  diffusion  layer  (GDL)  [2,8-10].  Man¬ 
ufacturing  GDLs  to  have  mixed  hydrophobic/hydrophilic  pores,  in 
theory,  promotes  both  gas  transport  to  and  liquid  removal  from  the 
catalyst  layer  [2].  This  mixture  of  hydrophobic/hydrophilic  pores  is 
accomplished  by  treating  the  moderately  wetting  fibrous  carbon 
substrates  with  polytetraflourethylene  (PTFE  or  otherwise  known 
as  Teflon®),  where  an  empirical  rule  of  thumb  is  to  apply  5-30  wt% 
PTFE  for  optimal  water  management  in  GDLs  [9,11,12]. 
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One  of  the  key  properties  of  a  GDL  that  gets  modified  by  Teflon® 
treatment  is  the  capillary  pressure  versus  liquid  saturation  behavior 
[8].  Capillary  pressure  (Pc)  is  defined  here  as 

Pc  =  Pl~Pg  (1) 


where  PL  and  PG  are  the  liquid  and  gas  phase  pressures,  respectively. 
Capillary  pressure  is  a  thermodynamic  property  that  relates  the 
liquid-vapor  interface  to  surface  properties  and  pore  geometries 
via  the  modified  Young-Laplace  equation  [13]: 


Pc  =  - 


2 y  cos  0 

^Eff 


(2) 


where  y  is  the  surface  tension  of  the  liquid-vapor  interface,  0  is 
the  contact  angle  of  the  liquid  with  the  solid  (measured  within  the 
liquid  phase),  and  REff  is  the  effective  pore  radius  in  an  equivalent 
cylindrical  capillary.  For  liquid  water,  0  >  90°  defines  a  hydrophobic 
material  whereas  0  <  90°  is  a  hydrophilic  material.  From  Eqs.  (1 )  and 
(2),  one  sees  that  Teflon®  rich  regions  in  the  GDL  are  hydrophobic 
{0  ranging  from  98°  to  112°  for  water  on  PTFE  [7,14,15])  and  are 
expected  to  have  a  positive  Pc,  whereas  the  carbon  rich  regions  are 
hydrophilic  (with  0  reported  to  be  80°  [7]  or  86°  [14]  for  water  on 
graphite)  and  are  expected  to  display  a  negative  Pc. 
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GDL  capillary  pressure  measurements  at  different  liquid  sat¬ 
urations  (Sl)  have  shown  a  mixture  of  both  hydrophilic  and 
hydrophobic  surfaces,  even  when  no  Teflon  is  added  to  the  structure 
[16-22],  where  liquid  saturation  is  defined  as 

c  J4_  =  _A_  (3) 

1  Wore  ^Wotal  1  ^ 

with  Vi  the  liquid  volume,  Wore  the  pore  volume,  e  the  porosity, 
and  Wotai  the  total  material  volume.  The  resulting  GDL  capillary 
pressure  vs.  liquid  saturation  curves  (Pc  vs.  Sl,  or  Pc(Sl))  provide  a 
constitutive  relationship  describing  the  integrated  capillary  prop¬ 
erties  (internal  wetting  and  pore  geometry)  in  the  porous  material. 
Additionally,  these  Pc(SL)  curves  are  dependent  on  the  history  and 
whether  the  liquid-vapor  interface  is  advancing  (liquid  intrusion) 
over  an  unwetted  surface  or  receding  (gas  intrusion)  over  a  wetted 
surface,  otherwise  known  as  hysteresis  [13,23-25]. 

Pc(Sl)  measurement  is  an  emerging  area  of  research  in  GDLs. 
Gostick  et  al.  [16]  used  contact  porosimetry  [26]  and  measured  the 
PC(SL)  curves  for  gas  intrusion  using  water  in  hydrophilic  pores  only. 
Kumbur  et  al.  [17-19]  also  used  contact  porosimetry  to  measure 
PC(SL)  curves  for  gas  intrusion  using  water  in  hydrophilic  pores  but 
also  looked  at  the  effects  of  compression  and  temperature.  Fair- 
weather  et  al.  [20]  developed  a  microfluidic  device  and  was  the 
first  to  measure  the  Pc(Sl)  curves  for  both  liquid  and  gas  intru¬ 
sion  while  probing  both  hydrophilic  and  hydrophobic  pores.  Gostick 
et  al.  [21]  recently  used  a  similar  method  to  Fairweather  et  al. 
to  probe  both  hydrophilic  and  hydrophobic  pores  by  controlling 
the  gas  phase  pressure;  they  observed  the  same  hysteresis  when 
cycling  between  liquid  and  gas  intrusion.  Nguyen  et  al.  [22]  used  a 
volume-displacement  technique  and  optically  tracked  the  change 
in  liquid  volume  but  did  not  observe  any  hysteresis  in  their  PC(SL) 
curves.  Although  Pc(Sl)  curves  provide  a  means  of  determining  the 
hydrophilic  and  hydrophobic  character  of  the  material,  the  specific 
pore  size  and  internal  wetting  properties  are  convoluted  through 
Eq.(2). 

Weber  et  al.  [7]  is  the  first  to  combine  independent  pore 
size  distribution  data  with  the  mixed  wetting  properties  of  the 
different  materials  present  in  GDLs.  Their  model  consisted  of  a 
bundle  of  capillaries  with  the  GDL  pore  size  distribution  deter¬ 
mined  from  mercury  intrusion  porosimetry  (MIP)  [27,28]  data. 
Each  pore  was  randomly  assigned  a  single  discrete  hydrophilic 
or  hydrophobic  contact  angle  in  proportion  to  the  amount  of 
hydrophilic/hydrophobic  material  present  in  the  GDL.  This  model 
was  then  capable  of  determining  the  extent  of  saturation  at  differ¬ 
ent  capillary  pressures,  with  PC(SL)  curves  showing  both  wetting 
and  non-wetting  characteristics.  Despite  the  significant  advance 
represented  by  Weber’s  work,  the  use  of  single-valued  contact 
angles  representing  the  wetting  and  non-wetting  surfaces  of  the 
material  produced  an  aphysical  feature,  namely  a  discontinuity 
at  liquid  saturations  where  all  the  hydrophilic  pores  were  filled. 
Recently,  Sinha  et  al.  [6]  used  a  uniform  distribution  of  contact 
angles  to  understand  the  behavior  of  water  in  GDLs  having  mixed 
wetting  properties.  A  uniform  contact  angle  distribution  is  a  com¬ 
putationally  straightforward  way  to  remove  aphysical  features  in 
the  capillary  behavior  of  GDLs,  but  the  most  realistic  representation 
is  likely  to  be  somewhere  between  Weber’s  [7]  single  values  and 
Sinha’s  [6]  uniform  values.  Specifically,  a  material  with  chemical 
heterogeneous  surfaces  will  show  an  averaged  wetting  described 
by  the  Cassie-Baxter  equation  [14,29,30]: 

COS  6>0bs  =/i  COS  0j  +  (1  -  ft )  COS  02  (4) 

where  0Obs  is  the  observed  contact  angle,  6\  and  02  the  contact 
angle  of  the  two  different  surfaces,  and  /i  the  fractional  coverage 
of  surface  1.  Moreover,  intrinsic  roughness  adds  variability  to  the 


Table  1 

Physical  properties  for  the  24  mm  diameter  Toray  TGP-H-090  GDL  sample. 


Property 

Value 

Thickness  (L)  (pm) 

280a 

Porosity  ( e ) 

0.78a 

Wore  (pL) 

99 

Advancing  contact  angle,  external  (#Adv,Ext)  (°) 

153° 

Receding  contact  angle,  external  (PRec,Ext)  (°) 

47° 

a  Denotes  manufactured  supplied  properties. 


contact  angle,  as  described  by  Wenzel’s  equation  [14,24,29,31  ]: 

cos  P0bs  =  r  cos  0o  (5) 

where  r  is  the  rugosity  or  roughness  factor  and  0o  is  the  intrin¬ 
sic  contact  angle.  Thus,  for  a  GDL  with  chemical  heterogeneity  and 
roughness  generated  by  Teflon®  or  binder  [11]  that  is  randomly 
distributed  throughout  the  porous  geometry,  one  would  expect  a 
Gaussian  distribution  of  contact  angles. 

In  this  work  we  present  a  method  for  analyzing  internal  surface 
wetting  in  GDLs  using  independent  data  from  MIP  and  PC(SL)  mea¬ 
surements  via  a  bundle  of  capillaries  model  that  assumes  a  realistic 
Gaussian  wetting  distribution.  This  approach  carefully  deconvo- 
lutes,  for  the  first  time,  the  intrinsic  material  properties  from  the 
porous  geometry  in  a  manufactured  GDL,  providing  the  kind  of  data 
that  can  guide  better  GDL  designs. 

2.  Experimental  methods 

2.1.  Materials 

A  commercially  available  as-received  untreated  Toray  TGP-H- 
090  GDL  sample  (Toray  CFA,  Flower  Mound,  TX)  was  investigated. 
External  advancing  and  receding  contact  angle  measurements  were 
performed  with  a  Rame-Hart  Model  100  contact  angle  goniome¬ 
ter  to  determine  the  extent  of  hysteresis.  Table  1  lists  the  physical 
properties  of  the  sample. 

2.2.  Microfluidic  system 

A  schematic  of  the  apparatus  used  to  measure  capillary  pressure 
in  the  GDL  samples  is  shown  in  Fig.  1.  A  liquid  manifold  (labeled 

A)  connects  a  syringe  pump  to  the  GDL  sample  housing  (labeled 

B)  and  a  differential  pressure  sensor  (labeled  Pc).  A  gas  manifold 
(labeled  C)  connects  the  top  of  the  sample  housing  to  a  second  pres¬ 
sure  sensor  (labeled  Pq)  and  two  gas  lines  providing  humidified  air 
and  a  vacuum  source.  The  faces  between  the  two  manifolds  and 


Fig.  1.  The  cell  used  to  measure  capillary  pressure  in  GDL  materials.  A  syringe  pump 
is  used  to  control  liquid  movement  into  and  out  of  the  GDL  sample.  A  liquid  manifold 
(A)  connects  the  pump  to  the  bottom  of  the  GDL  sample  housing  (B),  and  a  gas 
manifold  (C)  connects  the  top  of  the  sample  to  lines  for  vacuum  and  humidified  air, 
selected  by  valve  VI.  A  pressure  sensor  (Pc)  measures  the  difference  between  the 
liquid  and  gas  pressures.  A  second  sensor  (PG)  and  a  needle  valve  (V2)  are  used  to  set 
the  vacuum  pressure  during  the  initial  fill  procedure.  A  heating/cooling  recirculator 
passes  temperature  control  fluid  through  the  liquid  manifold. 
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Fig.  2.  Pore  size  distribution  of  Toray  TGP-H-090  showing  the  incremental  volume 
(AV)  contained  at  each  effective  radii  (PEff,i)- 


the  sample  housing  are  sealed  by  o-rings.  The  GDL  sample  itself  is 
sandwiched  between  a  hydrophilic  membrane  on  the  liquid  side 
and  a  hydrophobic  membrane  on  the  gas  side  to  provide  a  gas  bar¬ 
rier  and  liquid  barrier,  respectively.  The  construction  details  of  the 
GDL  sample  housing  have  been  described  previously  [20]. 

2.3.  Operation 

Before  each  experiment,  the  cell  is  filled  under  vacuum;  the  aspi¬ 
rator  line  is  selected  with  a  three  way  valve  (labeled  VI  in  Fig.  1 ),  the 
air  pressure  is  set  to  ~4  kPa  absolute  using  the  needle  valve  (labeled 
V2),  and  the  liquid  side  of  the  cell  is  filled  with  degassed  water  from 
the  syringe  pump.  The  liquid  manifold,  hydrophilic  membrane,  and 
GDL  sample  are  all  wetted,  while  the  hydrophobic  membrane  pro¬ 
vides  an  effective  barrier  to  the  water.  When  the  GDL  sample  is 
saturated  indicated  by  Pc  >  30  kPa,  the  vacuum  is  released,  and  a 
testing  procedure  is  followed  similar  to  one  previously  reported 
[20].  To  evenly  wet  the  GDL  sample,  the  syringe  pump  repeatedly 
cycles  liquid  water  into  and  out  of  the  GDL  sample  at  a  constant 
flow  rate  of  10  pi  min-1 .  During  cycling  the  difference  between  the 
liquid  and  gas  pressures  is  measured  using  the  “Pc”  pressure  sensor, 
with  the  volume  limits  chosen  such  that  the  measured  Pc  oscillates 
between  -30  and  +30kPa.  After  at  least  three  constant  rate  cycles, 
a  stepwise  cycle  is  performed  and  the  Pc(Sl)  curve  is  determined  as 
described  previously  [20].  Air  is  humidified  by  passing  it  through 
a  gas  sparge  bottle  and  supplied  to  the  gas  manifold,  to  minimize 
evaporation  from  the  sample.  Cell  temperature  is  controlled  by  a 
heating/cooling  recirculator  bath  attached  to  a  heat  exchanger  built 
into  the  liquid  manifold.  All  the  data  included  here  were  collected 
at  25  °C  for  three  repeated  measurements. 

2.4.  Mercury  intrusion  porosimetry  (MIP) 

Mercury  intrusion  porosimetry  (MIP)  was  performed  on  the 
GDL  sample  using  a  Micromeritics  AutoPore  IV  9500  (Micromeritics 
Instrument  Corporation,  Norcross,  GA).  The  pressure  range  applied 
ranged  from  0.4  to  30,000  psia.  At  each  applied  pressure,  the  incre¬ 
mentally  intruded  volume  ( AV)  of  mercury  was  recorded. 

3.  Results  and  analysis 

3.1  Pore  size  distribution  from  MIP  data 

Fig.  2  shows  the  incremental  volume  ( AV)  associated  with  each 
effective  pore  radius  {REff,i)  acquired  from  MIP  data  for  the  Toray 
TGP-FI-090.  Following  standard  procedures,  a  bundle  of  capillaries 


model  was  used  to  convert  each  applied  pressure  to  an  REfft/  via 
Eq.  (2)  by  assuming  6  =  130°  [27].  The  majority  of  the  GDL  volume 
occurs  in  pores  of  radii  near  10  p,m.  Dividing  each  A  V  data  point  in 
Fig.  2  by  the  individual  pore  volume  (tvR^.L)  at  that  effective  pore 
radius  provides  an  estimate  of  the  number  of  capillaries  in  that 
size  range.  Such  calculations  reveal  that  large  pores  (ftEff ti  ~  10  |xm) 
occupy  most  of  the  volume  but  are  a  lot  less  prevalent  in  num¬ 
ber  than  smaller  pores  (ftEff  *  ~  0.01  pan),  despite  the  fact  that  small 
pores  occupy  very  little  volume.  We  use  the  data  in  Fig.  2  directly 
in  our  bundle  of  capillaries  model  to  defines  the  RE fffl-  found  in  our 
materials.  Thus,  by  attaching  realistic  wetting  properties  (i.e.,  con¬ 
tact  angles)  to  these  capillaries,  a  Pc(Sl)  curve  can  be  generated. 

3.2.  Contact  angle  distribution  and  the  Pc(Sl)  curve 


The  contact  angle  distribution  is  described  by  a  normalized 
Gaussian  distribution  that  relates  the  probability  (xjr)  of  finding  a 
pore  with  a  contact  angle  0j, 

.  ,a  ,  exp{-(0j  -0Mean)2/2cr2} 

Wj)  =  ' - 5 -  (6) 

E^expH^  -  #Mean)  /2<T2) 

where  0Mean  is  the  mean  contact  angle,  a  is  the  standard  deviation  of 
the  contact  angle  due  to  heterogeneity  and  roughness,  and  M  is  the 
total  number  of  discrete  contact  angles  selected  for  computational 
purposes.  For  all  our  calculations,  M=  181  corresponding  to  0j  rang¬ 
ing  from  0°  to  180°  separated  by  integer  values.  We  further  assume 
that  %(,,  and  0j  are  uncorrelated,  so  that  we  can  randomly  assign 
contact  angles  among  all  pores  computed  through  independently- 
acquired  MIP  data. 

With  Eqs.  (2),  (6),  and  an  experimentally  measured  pore  size 
distribution,  we  are  now  in  a  position  to  determine  the  Pc  and  liquid 
saturation  level  for  every  possible  combination  of  PEff ?i  and  0j.  For 
any  given  capillary  pressure  we  select  {Pc  k ),  we  identify  a  filled  pore 
as  one  that  satisfies 


Pc(%f,i>  Qj)  <  PC,k 


(7) 


By  marching  through  every  possible  PEffl  and  0j  and  cataloging 
every  pore  satisfying  Eq.  (7),  the  saturation  of  the  GDL  at  the  given 
Pck  is  determined.  Algorithmically,  the  saturation  (SL)  is  deter¬ 
mined  by 


SdPC,k )  = 


(F  =  1  if  Pc(REff.i,£>j)sPc,k 
\  P  =  0  if  Pc(REff.i,  0j)  >  Pck 


(8) 


where  N  is  the  total  number  of  pores  and  F  is  unity  when  Eq.  (7)  is 
satisfied  for  a  pore  and  zero  otherwise.  By  systematically  increasing 
the  PCj<  and  determining  the  corresponding  SL,  a  PC(SL)  curve  is 
generated. 

One  way  to  use  the  capillary  bundle  model  is  to  specify  a  x/s  and 
determine  the  Pc(Sl)  curve.  Fig.  3 A  shows  two  sets  of  xj/  distribu¬ 
tions,  one  set  (labeled  i  and  ii)  has  a  constant  o  of  15°  but  two 
different  values  of  0Mean  (70°  and  110°)  and  the  other  (labeled  iii 
and  iv)  has  a  constant  0Mean  of  90°  but  two  different  values  of  a 
(15°  and  30°).  Fig.  3B  shows  the  PC(SL)  curves  that  result  after  using 
these  contact  angle  xfr  distributions  and  the  experimental  pore  size 
distribution  for  Toray  TGP-FI-090  given  in  Fig.  2.  The  Pc(Sl)  curve 
for  having  0Mean  =  70°  and  a  =  15°  (i)  shows  approximately  90%  of 
the  curve  being  hydrophilic,  consistent  with  0  =  90°  being  1.3  stan¬ 
dard  deviations  from  0Mean  =  7O°.  Similarly,  the  PC(SL)  curve  with 
^Mean  =  H0°  and  cr  =  15°  (ii)  shows  approximately  90%  of  the  curve 
being  hydrophobic.  Curves  (i)  and  (ii)  show  that  0Mean  mainly  shifts 
the  Pc(Sl)  curve  towards  either  higher  or  lower  Pc  values,  for  a  given 
SL.  High  capillary  pressure  regions  in  each  curve  correspond  to  the 
volume  of  the  sample  associated  with  small  capillary  radii. 

When  x/f  has  0Mean  =  90°,  curves  (iii)  and  (iv),  the  GDL  has  an  even 
split  of  hydrophilic  and  hydrophobic  pore  volume,  as  expected  for 
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Fig.  3.  (A)  Four  distributions  with  different  combinations  of  #Mean  and  cr: 
^Mean  =  70°,  (7  =  15°  (i),  0Mean  =  110°,  CT  =  15°  (ii),  $Mean  =  90°,  (7  =  15°  ( iii ),  and  $Mean  =  90°, 
(7  =  30°  (iv).  (B)  The  resulting  Pc(Sl)  curves  using  the  pore  size  distribution  of  Toray 
TGP-H-090  and  the  different  distributions  from  A. 


Contact  angle  (0)  /  deg 

Fig.  4.  (A)  The  experimental  Pc(Se,l)  ^ata  (symbols)  and  model  determined  best  fit 
Pc(Se,l)  curves  (lines)  for  Toray  TGP-H-090.  (B)The  resulting  distributions  for  both 
liquid  (solid  line)  and  gas  intrusion  (dashed  line). 


a  symmetric  contact  angle  distribution.  This  is  clearly  seen  in  the 
Pc(5 L)  curves  by  the  fact  that  Pc  =  0  kPa  occurs  at  SL  =  50%.  When  a  is 
increased  from  15°  (iii)  to  30°  (iv),  one  observes  that  cr  largely  con¬ 
trols  the  slope  of  the  PC(SL)  curve  in  the  middle  saturation  range. 
The  standard  deviation  in  the  contact  angle  distribution  ( a )  has  a 
large  impact  on  the  portion  of  the  sample  pore  volume  that  can 
generate  appreciable  capillary  pressure  deviations  from  the  mean 
capillary  pressure.  Specifically,  when  a  is  small  and  the  distribution 
of  contact  angels  is  narrow,  only  capillaries  with  small  PEff  can  gen¬ 
erate  appreciable  pressure  deviations,  but  as  cr  increases  and  the 
distribution  includes  a  wider  range  of  contact  angles,  the  fraction 
of  the  pore  volume  that  can  generate  high  capillary  pressures  also 
increases.  In  short,  the  model  results  in  Fig.  3  show  that  the  value  of 
^Mean  sets  the  overall  hydrophilicity/phobicity  of  the  GDL,  whereas 
cr  sets  the  fraction  of  the  pore  volume  that  deviates  significantly 
from  that  overall  character. 

3.3.  Best  fit  to  experimental  Pc(Se,l)  curves 

Instead  of  dictating  a  specific  ^  function,  the  most  representa¬ 
tive  can  be  determined  from  a  best  fit  of  the  capillary  model  to 
experimental  Pc(Sl)  data  using  0Mean  and  cr  as  parameters  of  the 
fit.  The  best  fit  was  determined  to  be  the  that  yielded  the  lowest 
root  sum  of  squared  error  (RSSE)  between  the  experimental  PC(SL) 
data  and  model.  The  best  fitting  protocol  was  performed  using  the 
Optimization  Toolbox  v.3.1  with  Matlab  R2007b.  The  details  of  the 
fitting  procedure  are  described  in  Appendices  A  and  B.  The  satu¬ 
rations  (Sl)  have  been  scaled  to  an  effective  saturation  (Se,l)  that 
ranges  from  0  to  1  in  order  to  remove  the  effects  of  the  residual  gas 


saturation  (SR)  and  irreducible  liquid  saturation  (Sir)  [13]: 


%l  = 


Sl-Sr 
Sir  -  Sr 


Sl  ~  SE|Min 

SL,Max  —  SE)Min 


(9) 


where  SL  Min  and  S^Max  are  the  minimum  and  maximum  liquid  sat¬ 
urations  that  correspond  to  the  SR  and  Sir,  respectively. 

Fig.  4A  shows  triplicate  experimental  Pc(SEil)  data  points 
acquired  using  the  new  microfluidic  system  and  the  best  fit  Pc(SEil) 
curves  for  Toray  TGP-H-090  for  both  liquid  and  gas  intrusion.  Table  2 
presents  the  tabulated  values  of  0Mean  and  cr  from  the  best  fits  with 
the  corresponding  coefficient  of  determination  (R2 ).  The  calculation 
of  R 2  is  given  in  Appendix  A.  As  the  R 2  values  show,  both  the  liquid 
and  gas  intrusion  curves  fit  the  experimental  data  quite  well,  with 
81%  of  the  liquid  intrusions  experimental  variance  captured  by  the 
capillary  model  and  93%  of  the  gas  intrusion  data  variance  captured 
by  the  model.  Fig.  4B  shows  the  resulting  distributions  that  result 
from  the  best  fit  Pc(%l)  curves.  Because  the  liquid  intrusion  Pc(SEi l) 
data  shows  a  significant  volume  of  pores  being  both  hydrophilic  and 
hydrophobic,  the  resulting  distributions  are  expected  to  sample 
a  range  of  both  hydrophilic  and  hydrophobic  contact  angles.  Fur¬ 
thermore,  since  the  majority  of  the  pore  volume  is  associated  with 
pores  on  the  order  of  REfftl-  ~  10  [xm  (see  Fig.  2)  the  resulting  Pc(SEi E) 
curve  is  largely  attributed  to  these  pores  with  contact  angle  values 


Table  2 

Best-fit  parameters  for  Toray  TGP-H-090. 


$Mean  (°) 

o(°) 

R2 

Liquid  intrusion 

92 

10 

0.81 

Gas  intrusion 

52 

8 

0.93 
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described  by  the  *P  function  in  Fig.  4B.  Alternatively,  the  gas  intru¬ 
sion  Pc(%l)  data  show  a  majority  of  the  pores  being  hydrophilic. 
As  a  result,  is  centered  at  a  hydrophilic  contact  angle  with  low 
probability  of  finding  a  pore  with  0  >  90°.  Examining  the  gas  intru¬ 
sion  Pc(%l)  curve  further  shows  an  inflection  at  SEL~95%.  This 
occurs  because  of  the  discrete  combinations  of  PEff ti  and  0  that  can 
have  Pc  ~  0  kPa.  Based  on  Eq.  (2),  Pc  on  the  order  of  a  few  kPa  can 
only  be  achieved  by  having  either  large  PEff  l-  or  0~9O°.  Because 
the  distribution  for  gas  intrusion  does  not  extend  to  90°,  the 
pores  corresponding  to  Pc  ~  0  kPa  can  only  come  from  the  larger 
pores  (PEff  i  >  10  |jim)  (see  Fig.  2).  Thus,  the  majority  of  the  Pc(SEi E) 
curve  is  attributed  from  PEff ti  ~  10  pum  but  pores  with  PEff ti  ~  100  pum 
dominate  the  curve  at  high  saturations  (SEtL  >  90%). 

From  the  data  in  Fig.  4A,  we  clearly  observe  hysteretic  wet¬ 
ting  between  liquid  and  gas  intrusion  into  the  GDL.  Flysteresis  in 
the  Pc(Sl)  curves  was  also  observed  in  prior,  less  refined,  exper¬ 
iments  from  our  lab  [20],  and  by  Gostick  et  al.  [21].  Because  the 
pore  structure  of  the  GDL  material  is  unlikely  to  change  during  liq¬ 
uid  or  gas  intrusion,  we  assume  all  of  the  experimentally  observed 
hysteresis  in  the  Pc(Sl)  curves  comes  from  hysteresis  in  the  dis¬ 
tribution,  as  shown  in  Fig.  4B  and  Table  2.  Contact  angle  hysteresis 
is  observed  in  most  solid  wetting  processes  unless  extraordinary 
lengths  are  taken  in  surface  preparation  (smoothness,  chemically 
homogeneous),  regardless  whether  involving  internal  or  external 
surfaces  [14,15,24,29].  In  fact,  our  external  measurements  of  the 
advancing  and  receding  contact  angles  on  this  material,  Table  1 ,  also 
shows  strong  hysteresis,  as  do  prior  external  contact  angle  studies 
on  other  GDL  materials  [  11  ].  Interestingly,  hysteresis  arises  from  the 
same  chemical  heterogeneity  and  surface  roughness  that  we  use 
to  justify  our  Gaussian  distribution  of  wetting  angles  [14,15,24,29]. 
Thus,  the  hysteresis  in  capillary  properties  that  we  observe  here, 
previously  [20],  and  in  other  labs  [21]  is  further  evidence  of  the 
complexity  of  GDL  wetting. 

Bundle  of  capillary  models  are  known  to  have  certain  limita¬ 
tions.  For  example,  interpretation  of  MIP  data  is  affected  by  details 
of  the  porous  geometry,  such  as  pore  bodies,  throats,  and  connec¬ 
tions  [21  ],  and  these  features  are  not  fully  captured  in  simple  bundle 
of  capillary  models  [14].  Nonetheless,  bundle  of  capillary  models  are 
the  dominant  approach  for  interpreting  MIP  data  and  converting  it 
into  pore  size  distributions  because  the  model  captures  many  first- 
order  effects.  We  anticipate  similar  issues  from  our  use  of  bundle  of 


distribution,  we  have  shown  how  different  *P  distributions  can  give 
rise  to  different  Pc(Sl)  curves.  In  general,  the  mean  contact  angle 
#Mean  shifts  the  magnitude  of  the  Pc  values  and  the  contact  angle 
standard  deviation  o  controls  the  slope  and  position  of  the  ends 
of  the  Pc(Sl)  curves.  By  allowing  0Mean  and  o  as  free  parameters, 
we  found  the  representative  *P  distributions  that  allowed  the  bun¬ 
dle  of  capillaries  model  to  best  fit  experimental  Pc(%l)  data  for 
both  liquid  and  gas  intrusion  curve.  Additionally,  the  *P  distribu¬ 
tions  fall  reasonably  within  the  observed  external  contact  angles 
measured.  The  combined  experimental  and  modeling  analysis  we 
outline  here  provides  a  new  way  to  link  GDL  wetting  properties 
to  the  manufacturing  of  GDLs,  and  ultimately,  should  help  design 
membrane  electrode  assemblies  that  are  less  prone  to  flooding.  For 
example,  this  method  may  help  understand  the  fundamental  effects 
of  Teflon®  treatment  on  the  resulting  pore  structure  and  surface 
properties  in  order  to  better  determine  the  optimal  Teflon®  loading 
based  on  the  desired  properties.  To  galvanize  the  linkage  between 
GDL  processing-properties-performance,  we  have  begun  investi¬ 
gating  the  role  of  PTFE  treatment  and  the  effects  of  material  aging 
on  GDL  capillary  properties. 
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Appendix  A.  Scaled  sum  of  squares  error  and  overall 
goodness  of  fit 

A  weighting  is  attached  to  a  range  of  data  points  for  the  RSSE 
calculation.  This  is  performed  because  the  Pc  data  measured  at  Se>l 
around  0%  and  100%  have  lower  certainty  due  to  complications 
with  gas  compression  and  expansion,  respectively,  present  in  the 
microfluidic  cell.  As  a  result,  the  Pc  values  corresponding  to  these 
saturation  values  have  the  least  confidence.  Thus,  the  data  points 
within  a  range  A P  of  Pc(%l  =  50%)  are  weighted  by  a  factor  of  10 
while  the  remaining  data  points  are  unweighted  shown  by  Eq.  ( A.l ) 


RSSE  = 


£iox( 

k= 1 

v  Pc’k  1  Expt  ^-k  |  ] 

Model)  2  forPc  (50%)  -  AP  <  Pc  <  Pc  (50%)  +  A P 

t( ^ 

k=  1 

1  -  Pc  k  1 

1  Expt  1  Model 

^  for  all  else 

(A.1) 


capillary  models  to  interpret  capillary  pressure  data.  The  first-order 
effects  this  analysis  reveals  should  significantly  advance  our  ability 
to  diagnose  the  processing-structure-property  relationships  that 
are  essential  for  GDL  design  and  manufacturing. 

4.  Conclusions  and  implications 

We  have  shown  an  approach  to  deconvolute  the  influence  of 
pore  structure  and  internal  wetting  properties  from  Pc(SEt L)  mea¬ 
surements  using  a  bundle  of  capillaries  model.  This  work  is  the 
first  to  report  representative  contact  angle  (*P)  distributions  that 
describe  the  internal  wetting  properties  of  Toray  TGP-FI-090  for 
both  liquid  and  gas  intrusion  that  best  fit  experimental  PC(SE> L) 
data.  In  our  model  we  have  taken  into  account  the  chemical  het¬ 
erogeneity  and  surface  roughness  of  the  GDL  materials  by  assuming 
a  single  Gaussian  distribution  to  describe  *P.  For  a  given  pore  size 


with  L  being  the  total  number  of  experimental  Pc  data  points.  It  was 
found  that  A P  of  4  kPa  consistently  captured  the  data  with  the  most 
confidence  quite  well.  Because  the  scaling  procedure  above  weights 
certain  points  more  heavily  than  others,  based  on  experimental 
observations,  it  is  important  to  ensure  there  remains  an  overall 
high  quality  fit.  The  overall  goodness  of  these  fits  is  assessed  via 
an  unweighted  coefficient  of  determination  (R2),  which  describes 
the  fraction  of  the  total  data  variance  captured  by  the  best- fit  model. 
It  is  computed  as 


^G  k  I  Expt  k  I  Model  J 

R2  =  1  -  k=X  L -  (A.2) 

ZXPc,fc|Expt-pc) 

k= 1 
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where 

L 

^  ^  Irxpt) 

Pc  =  — - 1 -  (A.3) 


Appendix  B.  Initial  conditions  for  best  fitting 

Due  to  the  highly  nonlinear  relationship  between  0Mean  and  a 
to  the  generated  PC(SL)  curve,  the  line  searching  process  using  the 
optimization  toolbox  tends  to  be  locally  trapped  near  the  initial 
value  at  a  local  minima.  In  order  to  have  higher  confidence  the  local 
minima  found  is  the  close  to  the  global  minimum,  different  initial 
values  spanning  the  range  of  each  parameter  was  used.  A  total  of  275 
combinations  were  used  to  cover  the  spectrum  of  possible  initial 
parameter  values. 
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